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Abstract 

The volume operator plays a pivotal role for the quantum dynamics of Loop Quantum Gravity (LQG), 
both in the full theory and in truncated models adapted to cosmological situations coined Loop Quantum 
Cosmology (LQG). It is therefore crucial to check whether its semiclassical limit coincides with the classical 
volume operator plus quantum corrections. 

In the present article we investigate this question by generalizing and employing previously defined 
coherent states for LQG which derive from a cylindrically consistently defined complexifier operator which 
is the quantization of a known classical function. These coherent states are not normalizable due to the 
non separability of the LQG Hilbert space but they define uniquely define cut ~ off states depending on 
a finite graph. 

The result of our analysis is that the expectation value of the volume operator with respect to coherent 
states depending on a graph with only rt— valent verticies reproduces its classical value at the phase space 
point at which the coherent state is peaked only if n = 6. In other words, the semiclassical sector of LQG 
defined by those states is described by graphs with cubic topology! This has some bearing on current 
spin foam models which are all based on four valent boundary spin networks. 
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1 Introduction 



The volume operator for Loop Quantum Gravitjo (LQG) enters the quantum dynamics of the theory in a 
prominent way. Without it, Hamiltonian constraint operators I4j, Master constraint operators [5] or physical 
Hamiltonian operators [6j cannot be defined. This is also true for truncated models of LQG such as Loop 
Quantum Cosmology (LQG) [7] which is supposed to describe well the homogeneous sector of the theory. 
In LQG one generically finds that the quantum evolution of operators corresponding to classically singular 
Dirac observables remain finite. This feature can be traced back to the way that inverse powers of the local 
volume, which enter into the expression for the triad, cotriad and other types couplings between geometry 
and geometry or geometry and matter, are quantized in the full theory [3]. Namely such derived operators 
are obtained as commutators between holonomy operators and powers of the volume operator. 

In view of its importance, it is of outmost interest to verify that the classical limit of the LQG volume 
operator coincides with the classical volume. By this we mean that the expectation value of the volume 
operator with respect to suitable semiclassical states which are peaked on a given point in phase space 
coincides with the value of the classical volume at that phase space point up to small corrections and that 
its fluctuations are small. 

Notice that there are actually two versions of the volume operator discussed in the literature [8l l9] which 
result from nonequivalent regularisations of the products of operator valued distributions that appear at 
intermediate stages. However, only the operator in [9,j survives the consistency test of [10], namely that 
writing the volume in terms of triads which then are quantized using the above mentioned commutator, 
delivers the same operator up to h corrections as the direct quantization. This check is not unimportant as 
otherwise we should not trust the triad and cotriad quantisations that enter the quantum dynamics. 

The semiclassical analysis of the volume operator has not been carried out yet although in principle 
suitable semiclassical (even coherent) states for LQG are available [11]. The reason for this is that the 
spectral decomposition (projection valued measure) of the volume operator cannot be computed exactly in 
closed form which however is needed in order to do exact practical calculations. More precisely, the volume 
operator is the fourth root of a positive operator Q whose matrix elements can be computed in closed form 
|12j but which cannot be diagonalized analytically. More in detail, the volume operator has discrete (that is, 
pure point) spectrum and it attains a block diagonal form where the blocks are labelled by the graphs and 
spin quantum numbers (labelling the edges of the graph) of spin network functions (SNWF) [13j. SNWF's 
form a convenient basis of the LQG Hilbert space pUJ which is the unique (cyclic) Hilbert which carries 
a unitary representation of the spatial diffeomorphism group and the Poisson*— algebra of the elementary 
(flux and holonomy) variables [T5]. The blocks turn out to be finite dimensional matrices whose matrix 
elements can be expressed in terms of polynomials of 6j symbols which fortunately can be avoidecH using a 
telescopic summation technique [16] related to the Elliot ~ Biedenharn identity [T7j so that a manageable 
closed expression results. However, the size of these matrices grows exponentially with increasing spin 
quantum numbers and since the expression for coherent states is a coherent superposition of SNWF's with 
arbitrarily large spin quantum numbers, a numerical computation of the expectation value using numerical 
diagonalisation techniques which are currently being developed [18] is so far out of reachU 

The way to make progress is to use semiclassical perturbation theory developed in [19] and already 
applied in [20l[2T]. The idea is quite simple: In practical calculations one needs the expectation value of Q"^ 
where g is a rational number in the range < q < j- Introduce the "perturbation operator" X := — 1 
where the expectation value < Q > of the positive operator Q is exactly computable. Notice that X is 
bounded from below by —1. Then trivially < Q'' >=< Q < [1 + X]'^ >. Now we use that there exist 
positive numbers p such that the classical estimates 1 + qx — px^ < (1 + xY < 1 qx are valid for all 

^See [Il[2] for recent books and [3] for recent reviews. 

^Notice that Racah's formula provides a closed expression for the 6j symbol but that it involves implicit sums and factorials 
involving large integers which quickly becomes unmanageable even numerically. 

''The coherent superposition contains a damping factor which suppresses large spins and thus large matrices so that one 
may truncate the involved infinite series over spin quantum numbers at finite values making only negligible errors, but still the 
computational effort is currently too high for supercomputers, see e.g. the computation time estimates in reported [18j . 
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X > —1. Finer estimates of this form involving arbitrary powers of x are also available |20J. In view of the 
spectral theorem, this classical estimate survives at the quantum level and we have Y- < Y < 1+ where 

y+ = 1 + qx, y_ = y+ - px^, y = (i + x)'^. it follows that < y >e [< y+ > -p < x^ >, < y+ >]. 

However, < y+ >= 1 and < >= [< > — < Q >^]/ < Q >^ is proportional to the relative fluctuation 
of Q which of order of h It follows that to zeroth order in h we may replace < Q'' > by < Q which 
is computable and whose error estimates are also computable as shown above. 

While possible, the exact computation of < Q >, < >, ..is still quite involved. Here the following 
result, proved in [11], is convenient: The computation of these expectation values (more generally for low 
order polynomials in flux operators) for SU{2) spin network states coincide, to zeroth order in h, with the 
corresponding calculations for U{1)^ spin networks. On those, the volume operator is even diagonal. Hence 
we conclude that, as long as we are only interested in the zeroth order in h contribution, we may evaluate 
the expectation value of the volume operator for a fictive theory in which we may replace the non Abelian 
group SU{2) by the Abelian group (Ul)^ which makes all calculations dramatically simpler. 

The coherent states developed for LQG so far all have been constructed using the complexifier method 
reviewed in [22] which generalizes the coherent state construction for phase spaces that are cotangent bundles 
over compact groups developed in [23pl . This involves the heat kernel evolution of the delta distribution 
(which is the matrix element of the unit operator in the Schrodinger (position) representation) with respect 
to a generalized Laplace operator, called the complexifier, followed by a certain analytic continuation. Now 
the unit operator in LQG can be written as a resolution of unity in terms of SNWF's and although the heat 
kernel is a damping operator, since the SNWF are not countable (the LQG Hilbert space is not separable), 
the resulting expression is not normalizable. It defines a well defined distribution (in the algebraic dual of 
the finite linear span of SNWF's) which can be conveniently written as a sum over cut - off states labelled 
by finite graphs. These states, called "shadows" in [23], are normalizable and one can use those in order to 
perform semiclassical calculations. Of course, one expects that only cut - off states labelled by graphs which 
are sufficiently fine with respect to the classical three metric to be approximated are good semiclassical 
states. 

Hence we see that the input in the semiclassical calculations consists in the choice of a complexifier and 
the choice of a graph. One may wonder why the complexifier has to be chosen and is not dictated by the 
dynamics of the theory, that is, that the coherent states remain coherent under quantum time evolution 
(e.g. for the harmonic oscillator, the complexifier is strictly the Laplacian on the real line). The reason is 
twofold: 

On the one hand one could perform constraint quantization and then we are working at the level of the 
kinematical Hilbert space on which the quantum constraints have not been imposed yet and all we want to 
make sure is that the volume operator, which is not gauge invariant (and thus does not preserve the physical 
Hilbert space) has a good classical limit on the kinematical Hilbert space (we would not trust a constraint 
quantization for which not even that was true). 

On the other hand, one could also work at the level of the physical Hilbert space as outlined in [6] and 
then one would expect that the time evolution of any reasonable choice of complexifier coherent states with 
respect to the physical Hamiltonian defined there remains coherent (and peaked on the classical trajectory) 
for sufficiently short time intervals. Notice that for interacting theories as simple as the an-harmonic 
oscillator globally stable coherent states have so far not been found. 

The choice of the complexifier will be guided by practical considerations, namely that it be diagonal on 
SNWF's and that it is a damping operator which makes the heat kernel evolution of the delta distribution 
restricted to a graph normalizable. Moreover, it should be gauge invariant under the SU (2) gauge group. As 
far as the choice of the graph is concerned, for practical reasons one will choose graphs that are topologically 
regular, that is, have constant valence for each vertex. Indeed, the semiclassical calculations performed in 
[20l [21] were done using graphs of cubic topology with good semiclassical results. 

The existing literature on such coherent states for LQG can be divided into two classes, depending on a 
certain structure that defines them: 

''See also [24j for related ideas valid for Abelian gauge theories such as Maxwell theory or linearized gravity. 
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On the one hand, there are gauge covariant flux coherent states which depend on cohections of surfaces 
and path systems inside them [H]. On the other hand, there are area coherent states which only depend 
on collections of surfaces but involve area operators rather than flux operators [22]. We will review and 
generalize both constructions. In its most studied incarnations, the collection of surfaces involved in [11] is 
defined by a polyhedronal partition of the spatial manifold while the collection involved in [22] is in terms 
of a parquette of foliations of the spatial manifold. 
The novel results of this article are as follows: 

1. Cylindrically consistent complexifier 

The important thing to notice is that both constructions are cylindrically consistent: Both derive from 
a complexifier which is graph independently defined in terms of the surfaces (and path systems) only 
and has consistent cylindrical projections to spin network functions over any graph. We stress this, 
because from the discussion in ^llj one may get the impression that the states constructed there do 
not come from a single complexifier operator. We clarify this in this paper by explicitly constructing 
those projections. The apparent contradiction is resolved by observing that in [llj cut - off states 
were displayed only for graphs which are dual to a fixed polyhedronal complex. In this paper we treat 
the general case which is more complicated. We also give a compact expression for the eigenvalues of 
the complexifier of type [22] which is missing there. 

2. Graphs of arbitrary valence 

The real motivation for the present article is that the calculations performed in [20^ [2T] raised the 
suspicion that the volume operator attains acceptable coherent state expectation values only for graphs 
of cubic topology. Therefore in this paper we perform the semiclassical analysis with respect to 
topologically regular graphs of valence n with n = 4, 6, 8. Ideally one would like to treat the general 
case but this leads to complicated book keeping problems. Fortunately, the cases we consider turn out 
to be sufficient to answer the afore posed question. 

The outcome of our analysis is that, no matter whether one uses the states [H] or [22], the correct semiclas- 
sical limit is attained with these states for n = 6 only. In other words: 
Up to now, there are no semiclassical states known for other than cubic graph topology! 

We interpret this result as saying that the states we constructed are simply not semiclassical for the 
volume operator unless n = 6. The ratio consisting of expectation value of the volume operator with respect 
to the coherent states for graphs of valence n ^ 6 divided by the classical value deviates from unity by a 
geometrical factor g„ depending on n which is of order unity and not at all of order h. One might argue that 
therefore one should rescale the volume operator by 1/qn- However, first of all this is not allowed because 
the normalization is fixed [10] and secondly even if one would rescale, the statement would still be that the 
correct expectation value is attained for a unique valence only. 

It may be that there are other semiclassical states coming from a different complexifier or from an entirely 
different method which does not have this problem. However, notice that all coherent states that were 
constructed so far for free field theories on Minkowski space come from a complexifier, so the complexifier 
method is natural and well tested in other contexts [22j. Secondly, as shown in [22], the choice of a SU(2) 
invariant and background metric invariant complexifier for LQG with computable spectrum and the required 
damping behavior leaves relatively little possibilities and to the best of our knowledge, the states proposed 
in [11^ 122] and used here are the simplest and only ones known so far. Therefore, even if better suited states 
exist, they will be very hard to guess and even harder to do practical calculations with. 

The implication of our result for LQG then seems to be that the semiclassical sector of the theory is 
spanned by SNWF based on cubic graphs. This has some bearing for spin foam models [26J which are 
supposed to be - but have so far not been proved to be - the path integral formulation of LQG. Spin foams 
are certain state sum models based on simplicial triangulations of four manifolds whose dual graphs are 
therefore five valent. Since the intersection of this graph with a boundary three manifold is therefore four 
valent, we see that spin foam models based on simplicial triangulations corresponds to boundary Hilbert 
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spaces spanned by spin network states based on four valent graphs onl}0. 

But even if these mismatches between LQG and spin foams could be surmounted, the result of our 
analysis seems to say that the boundary Hilbert space of current spin foam models does not contain any 
semiclassical states! This seems to contradict recent findings that the graviton propagator derived from spin 
foam models comes out correctly [311. However, notice that these results only show that the propagator 
comes out with the correct fall off behavior while the correct tensorial structure has not been verified yet. 
An easy way to possibly repair this is to generalize spin foam models and allow for arbitrary, in particular 
cubic, triangulations as suggested in [32l [33]. 

The outline of this paper is as follows: 



In section two we review the general complexifier method and generalize the complexifier coherent states 
proposed in jlll I22j . In particular, we review the two families of coherent states introduced so far which are 
based on complexifiers constructed from squares of flux or area operators respectively. We also explicitly 
compute the corresponding cut - off states for SU(2) and show when and why it is justified to work with a 
Active U{1)^ theory instead as far as the zeroth order in h for the expectation values is concerned. 

In section three we construct explicitly convenient tetrahedronal, cubical and octaheronal cell decom- 
positions of the spatial manifold which also define regular dual n = 4, 6, 8 valent dual graphs. These 
constructions are needed for the explicit calculations in section four and for our companion paper |34j . 

In section four we compute the expectation value of the volume operator for the constant valence cases 
n = 4, 6,8 for the states of type [11]. We will do the same in our companion paper [34J for the states of 
type [22]. In addition, in [3l] we answer the following question: While the flux complexifier is built from a 
discrete set of fluxes, the area complexifier is built from a continuous set of areas whose underlying surfaces 
fill all of space. While this makes the calculations more cumbersome, it might improve the quality of the 
corresponding coherent states in the sense that the expectation values do not strongly depend on the cut 
- off graph chosen (as already mentioned above, at least when the graph has cubic topology). Hence, we 
examine the question whether the expectation value of geometrical operators is invariant under Euclidean 
transformations of the graph in case that the spatial metric to be approximated is Euclidean. 

In section five we summaries and conclude. 

2 Complexifier Coherent States 

In this section we review the complexifier method to construct coherent states. We will be brief, all details 
can be found in [22] and references therein. This section is divided into six parts: 

In the first we define the general complexifier framework. In the second we provide natural gauge invari- 
ant choices of complexifiers for background independent SU(2) gauge theories such as General Relativity 
and fictive U{1)'^ gauge theories. In the third we compute the associated coherent states. These are not 
normalizable because they involve a sum over an uncountable set of graphs but one can extract from them 
normalizable states in which the sum over all graphs is cut off to a finite number and which we display in the 
fourth part. In the fifth part we compare the properties of the states |11] and [22]. Finally, in the sixth we 
show that as far as the zeroth order in h is concerned, we may replace the complicated SU{2) calculations by 
the much simpler U{1)^ calculations under certain restrictions on the cut - off graph. This approximation 
is needed only for states of type [22]. For the states of type [11] it is sufficient to take as a cut - off graph 

^ As an aside, whether this boundary Hilbert space of spin foams really can be interpreted as the four valent sector of LQG 
is a subject of current debate even with the recent improvement [27] of the Barrett - Crane model [5S] mostly studied so far. 
There are two problems: First, the boundary connection predicted by spin foams does not coincide with the LQG connection 
[29] . Secondly, the four valent sector of the LQG Hilbert space is not a superselection sector for the holonomy flux algebra of 
LQG. In fact, the LQG representation is known not only to be a cyclic but even an irreducible representation [HO]- Therefore 
the four valent sector is not invariant under the LQG algebra. 
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any graph which is dual to the polyhedronal complex that defines the complexifier. 



2.1 General Complexifier Method 

We consider a symplectic manifold M = T*{C) which is a cotangent bundle over a configuration manifold 
C which may be infinite dimensional (we will suppress any indices in what follows). 

Definition 2.1. 

A complexifier C : M ^ is a sufficiently smooth, positive function on M with dimension of an action 
which has the following scaling behavior 

C{q,Xp) 

lim = oo (2.1) 

A^oo A 

where {q,p) are the canonically conjugate, real configuration and momentum coordinates on M. 
The reasons for these restrictions become evident in a moment. With the aid of C we define 

z := exp(-i{C, .})-q = Yl i^, q}{n) (2.2) 

n=0 

with inductively defined multiple Poisson brackets {C, /}(o) := /, {C,f}(^n+i) •= {C) {C*, /}(«)}• The 
meaning of "sufficiently smooth" is that all coefficients in the Taylor expansion (12. 2p exist. Notice that (12. 2|] 
defines a (complex valued) canonical transformation whence {z,z} = = (this is non trivial when 

dim(C) > 2). The scaling behavior implies that z, z can be used as coordinates for Ai. In fact, z E 
defines a complex polarization of Ai. 

We now assume that Ai can be quantized, that is, there is a representation {q,p) i— > {q,p) of the 
Poisson* algebra defined by {q,q} = {p,p} = 0, {p,q} = 1m] q = q, p = p on a Hilbert space of the form 
Ti. = L2{C, djj) where C = C in the finite dimensional case and otherwise C C C in the infinite dimensional case 
is a suitable distributional extension. That is, the operators satisfy (assuming careful domain definitions) 
[q^q] = [PjP] = 0, \p,q\ = ih !■}{] q^ = = P- Here C comes with some topology and // is a Borel 

measure on it. 

Assuming that also C has a quantization (7 as a positive, self adjoint operator (in field theories this is 
non trivial due to operator ordering and operator product expansion questions) we construct the operator 
representation of (j2.2p by substituting Poisson brackets by commutators divided by ih 

This formula explains the dimension restriction on C. The operator e^^^^ is well defined via the spectral 
theorem. We will refer to e~^/^ as the heat kernel and to z as the annihilation operator. 

The 5 distribution with support at q^ ^ C on the subset of Ti consisting of the continuous functions 
is defined by 

•^goM := V'(go) :=< <5go,V' >:= j_ dn{q) 5q^{q)^{q) (2.4) 

C 

where dq^ (q) is the integral kernel of the unit operator. We define for z G C the coherent state 

i,, := [e-^/%J,„_, (2.5) 

It is defined as "heat kernel evolution" followed by analytic continuation. In order that this makes sense, 
the function e~'^/'^5gp must not only be in Ji but also analytic in q^. This explains the positivity and 
scaling requirement on C which makes sure that the heat kernel is a damping operator such that at least 
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for separable TC the function (|2.5|) is not only normalizable but also analytic. Here we have assumed that 
the map C — > C*^; g i— > z in (j2.2p has an extension to some . 

The whole point of this construction is that one can easily verify 

z tjj^ = z tp;, (2.6) 

Thus, tl^z is an eigenfunction of the annihilation operators z which explains the notion "coherent state" . As 
is well known, property (j2.6p implies that the uncertainty relation for the self adjoint operators 

X := [z + z^]/2, y := -i[z - z^]/2 (2.7) 

is saturated, that is 

[< x2 >, -(< X >,)2] [< y2 >^ _(< y >^^)2] = 1 I < >J (2.8) 

where < . >z'=< ipz,- tpz > /W^PzW^ denotes the expectation value with respect to "02 (notice that tpz is in 
general not automatically normalized). This is a second property commonly attributed to coherent states 
[35]. 

Finally, under certain technical assumptions spelled out in [36j the completeness relation 

In = l_ dfi{qo) Sg, Sg,[.] (2.9) 
— c 

implies that there exists a measure on C such that 



ln= /„ diy{z)iPz <V'2,. > (2.10) 



c 



This concludes the general discussion. The interested reader may verify [22\ [M] that the coherent states for 
Maxwell Theory on Minkowski space result from the complexifier 

C=7T^/ d'^xdabE'^V^'^E^ (2.11) 
2k^ Jr3 

where E"' is the Maxwell electric field, A is the Laplacian on M^, k is the electric charge and a = Hk? is the 
Feinstruktur constant. 



2.2 Complexifiers for Background Independent Gauge Theories 

As explained in detail in [371 EH], gauge theories with compact gauge group G provide an almost perfect 
arena for the general theory summarized in the previous subsection. Let us explain this in some detail: 

1. Classical Phase Space 

The role of C is played by some space of smooth connections A over some £)— dimensional spatial 
manifold a. The role of M is then simply T*A. The configuration and momentum coordinates on 
this phase space then are simply real valued connection one forms Aa{x) (potentials) and Lie algebra 
valued vector densities Ej{x) (electric fields) respectively which enjoy the following Poisson brackets 

{Ai{x),AUy)} = {E^{x),EUy)} = 0, {Ef{x),AUy)} = ^ St 6^ S{x,y) (2.12) 

Here k denotes the square of the coupling constant of the gauge theory, a,b,c,.. = l,..,D denote 
spatial tensor indices and j, k, I, .. = 1, .., dim(G) denote Lie algebra indices. We will assume that G is 
connected, semisimple and take the convention that the internal metric is just 6jk- 
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2. Distributional Configuration Space 

Now consider arbitrary, finite piecewise analytic (more precisely semianalytic |15j) graphs embedded 
in a which we think of as collections of edges e, that is, piecewise analytic one dimensional paths which 
intersect at most in their endpoints, called the set V{'y) of verticies of 7. Denote by -^(7) the set of 
edges of 7. Given 7 consider functions cylindrical over 7 of the form 

/ : ^ ^ C; ^ ^ f{A) = mA{e)},^Ei-y)) (2.13) 

where : G'^WI ^ C is a complex valued function on |-E(7| copies of G and A{e) denotes the holon- 
omy of A along e. Functions of the form ()2.13p form an Abelian * algebra under pointwise operations 
with the involution given by complex conjugation. We can turn it into an Abelian C*— algebra, usually 
called Cyl (cylinder functions) with respect to the sup norm on A that is 

ll/lh=sup \f{A)\ (2.14) 

As is well known [39], Abelian C*— algebras 21 are isometric isomorphic to another Abelian C*— algebra 
which consists of continuous functions on a compact Hausdorff space A (21), called the spectrum of 21. 
Denote the spectrum of Cyl by A. It has a nice geometrical interpretation as a space of generalized 
connections in the sense that the holonomy of ^4 G ^ satisfies all the usual algebraic relations satisfied 
by smooth holonomies, that is A{eoe') = A{e)A{e') if the end point of e is the beginning point of e' and 
A{e~^) = (^(e))~^ , but that smoothness or even continuity is no longer required. The topology on A 
is the Gel'fand topology which in this case boils down to saying that a net of generalized connections 
converges when the corresponding net of holonomies for all possible paths converges. See [14^ HOt [2] 
for more details. 

3. Hilbert Space 

Being a compact Hausdorff space, a natural set of representations of the Poisson*— algebra generated 
by all the holonomies and all the electric fluxes through co-dimension 1 (piecewise analytic) surfaces 
S 

Ej{S) := [ *Ej (2.15) 
Js 

where *Ej denotes the pseudo {D — l)form dual to Ej, should be of the form 7i = L2{A,dn) where 
fi is a Borel probability measure. It turns out that all cyclic representations that carry a unitary 
representation of the diffeomorphism group Diff(o") are of this form |T5] and the corresponding measure 
is unique and was first discovered in [14]. See e.g. [2\ for all details. For our purposes it is enough 
to know that TC admits a natural orthonormal basis, called spin network functions (SNWF). They are 
labelled by a graph 7, a collection vr = {'n'e}e(^E{-y) of irreducible, non trivial representations of G and 
collections m = {me}eeE{-j)^ ^ = {^e}ee-E(7) of integers rUe, = 1, ..,dim(7re) labelling the edges of 
7 and are explicitly defined by 

T^,n,m,n{A) ■= H V^M^e) ^ (^(e) )]™.ne (2-16) 
ee£;{7) 

In order to further apply the general theory of the previous subsection we must provide a complexifier. The 
complexifier for Maxwell theory displayed in ()2.1ip is motivated by the fact that the associated annihilation 
operators are precisely those that enter the Maxwell Hamiltonian. In General Relativity there is no a 
priori Hamiltonian but there is the Hamiltonian constraint. Hence one might be tempted to choose a 
complexifier whose associated annihilation operator is related to the Hamiltonian constraint. Unfortunately, 
the Hamiltonian constraint is, in contrast to Maxwell theory, neither polynomial nor does it have a quadratic 
piece with respect to which a perturbation scheme can be defined. Hence, the notion of an annihilation 
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operator defined by the Hamiltonian constraint is ill defineci^. On the other hand, since we here just want 
to construct coherent states which approximate well our elementary holonomy and flux operators which are 
defined on the kinematical Hilbert space (on which the Hamiltonian constraint is not satisfied) in order 
to verify whether other kinematical operators (i.e. not invariant under the gauge motions generated by 
the spatial diffeomorphism and Hamiltonian constraints) such as the volume operator have been correctly 
quantized, the motivation to use the Hamiltonian constraint as a selection criterion for the complexifier is 
anyway less motivated. 

In lack of a better selection criterion, we take here a practical attitude and would like to consider a 
complexifier which comes close to the Maxwell one (j2.1ip which obviously satisfies all the requirements 
of definition 12.11 Since we have applications in General Relativity in mind, we must preserve background 
independence and therefore the Minkowski background Laplacian entering (|2.1ip must be replaced by some- 
thing background metric independent. One possibility is to use a background independent Laplacian which 
depends on the dynamical 3-metric of Qat which Ej det{E)\ is its triad. However, this would lead to 
a very complicated object with which no practical calculations are possible. In fact, the practical use of 
coherent states in Maxwell theory rests on the fact that (j2.1ip is quadratic in the momenta (electric fields) 
which leads to states which are basically Gaussians in both the position and the momentum representation. 
This motivates to keep our complexifier quadratic in the momenta as well. Furthermore, we must preserve 
G invariance. For Abelian gauge theories the electric fields are already gauge invariant but not for non 
Abelian gauge theories. 

Thus a first attempt would be to define as complexifier 

Coc / d''xqabE^Ei6^''/y^d^ (2.17) 

J a 

where had to replace the background metric b^h in (|2.1ip by the dynamical metric and in order to make 
()2.17p spatially diffeomorphism invariant we have included a density factor l/y^ det((7). However, it is easy 
to see that (|2.17p becomes 

C(xV= f d^x ^\det{E)\ (2.18) 

J a 

the volume functional. While it satisfies the requirements of a complexifier, and admits a quantization as a 
positive self adjoint operator, its spectral decomposition is not analytically available so that C = y is not 
practically useful. 

Hence, what we need is a gauge invariant, background independent expression, quadratic in the electric 
fields which preferably is non vanishing everywhere on a and which can be expressed in terms of (limits of) 
electric fiuxes since only those are well defined in the quantum theory. In [22] it is shown that in non Abelian 
gauge theories no quadratic complexifier based strictly on fluxes exists that meets all these requirements. 
The way out is to give up the the requirement that the complexifier is composed out of the fluxes but to 
allow more general objects than fluxes. There are basically two proposals in the literature. The flrst [11] 
replaces fluxes by gauge covariant fluxes. The second [22] replaces the flux by areas. We will review these 
two proposals separately. 

In what follows we assume that as in General Relativity the canonical dimension of E'^ is cm'' and 

that of J^a is cm~^ so that ()2.25p has dimension cm^^^. Since the kinetic term in the canonical action is 
Jjjdt J^d^EjAi/K it follows that Hk has dimension cm^^^. 

2.2.1 Gauge Covariant Flux Complexifiers 

Given a surface S select a point p{S) £ S. Furthermore, for each point x £ S choose a path ps{x) C S 
within S with beginning point p{S) and ending point x. Denote the path system by Vs- The gauge covariant 

^The situation slightly improves when a physical Hamiltonian is available, see [6]. 
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flux of E through S subordinate to the path system Vs and the edge es is defined by 

E,{S)t, := / AdAipsix))i{*E){x)) (2.19) 

Here irj are the Pauli matrices, *E = ^Cabcdx"" A dx^ ^'j^j ^'^'^ denotes the adjoint action of G on its Lie 
algebra. Obviously, (I2.19P transforms in the adjoint representation under gauge transformations at p{S). 

Let 5 be a collection of surfaces with associated path systems Vs for each S £ S and /x a measure on 
S. Let K be any positive definite, measurable function on S. A gauge covariant flux complexifier (GCFC) 
is defined by 

C := ^ dfiiS) K{S) [-^Tt{E{Sf)] (2.20) 

which is manifestly gauge invariant (one could absorb K, L into fi) . Here L is a parameter of dimension of 
length and we assume both fi, i^T to be dimensionless. 

The mostly studied case is when D = 3 and S = C^{V) is a discrete set of oriented surfaces which 
coincide with the faces (its sub 2-complex) of a polyhedronal cell partition P of a. In this case /x is just the 
counting measure and for convenience one chooses K = 1. We will denote the associated complexifier by 
Cp. In this case the complexified connection is given by 

Zi{x) = Aiix)-^ [ leabcdy'Ady' [Tr{E{S)Ad,^^,){T,))]6{x,y) (2.21) 

Notice that the series involved in Zi terminates at the first term. This is because when computing the 
second iterated Poisson bracket there is a double sum over surfaces involved but because the paths psix) 
are disjoint from 5' for S' ^ S there is no contribution from S" 7^ 5 to {Cp, Ai{x)}(^2)- S' = S there 
is in principle a contribution but by the regularization i2j the classical flux does not Poisson act on paths 
lying in its associated surface. 

This connection is distributional but fortunately we are only interested in the integral of (I2.2ip over one 
dimensional paths e given by 

iL^ [dx-[Zi{x)-Ai{x)]= E ^-('^'^)['^(^(^)Adps(-)(^^))] (2.22) 



S£C2{P) x&Sne 



where 



<Tz{S, e) = ^ ^ dx'' eabc dy'> A dy''5{x, y)5,,, (2.23) 

is the signed intersection number at z £ eCiS which here we have assumed to be an interior point (otherwise 
there is an additional factor of 1/2, in [2]). 

2.2.2 Area Compexifier 

Let 5 be a collection of surfaces, fj, a measure on S and K{S, S') a positive definite integral kernel. An area 
complexifier is given by the expression 

C = I diiiS) I dfi(S) K(S, S') Ar(S) Ar(S') (2.24) 

a 1^ Js Js 

where o is a parameter of dimension of length. Here Ar(5) is the gauge invariant "modulus of the electric 
flux" 



Ar{S) := / V^HFW) (2.25) 



which in General Relativity has the meaning of the area of S. 

The most studied case arises from a diagonal and constant integral kernel and the following choices of 
S and n respectively. 
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Definition 2.2. 

i. 

A stack s in a is a D- dimensional submanifold with the topology o/R x (0, 
ii. 

A stack family S = {sa} is a partition of a into stacks which are mutually disjoint. 
Hi. 

D families of foliations Fj, I = 1,..,D of a generated by vector fields d/dt^, I = 1,..,D are said to be 

linearly independent if the vector fields d/dt^ are everywhere linearly independent. 

iv. 

D stack families are said to he linearly independent provided that there exist D linearly independent foli- 
ations Fj such that the leaves of the foliation Fj is transverse to every stack in . That is, the intersection 

of any leaf Ljt, t G M of Fj with any stack in , called a plaquette, has topology (0, l]-^"-*-. 
V. The collection of the plaquettes s^ at is called a parquette at time t within Ljt. 

In general a will have to be partitioned into pieces that admit D linearly independent foliations each. 
We construct the complexifier for one such piece below, the complete complexifier is then the sum over the 
pieces. 

The complexifier defined by D linearly independent stack families is now is defined by 

Here we take the foliation parameter t to be dimensionless, a is a parameter with dimension cm^ so that 
C/his dimensionfree and = fl Ljt denotes the plaquette at time t within the stack in direction /. 
For Abelian gauge theories also the following simpler expression is available 

D 

I 

c 



^p^J2T.ldt[E,{pi,)]' (2.27) 



2Ka 

1=1 a 

which uses the gauge invariant flux rather than the areas. 

Let us now compute the complexified connections. Notice that due to the fact that each stack is foliated 
by squares with half open and half closed boundaries, for each x £ a and each direction / there exists a 
unique stack s^(a;) corresponding to a label ai{x) such that a; G s^. Likewise, for each direction / there exists 
a unique leaf Lit{x) corresponding to a time ti{x) such that x G Ln. Consider the one parameter family of 
embeddings X^^ : [0, 1)-^"^ — p^^, then there exists a unique ui{x) such that x = 
now set 

T ( \ — \A ( ^^at (^) \ I 

Jl[X) := |det( )\a=ai{x),t=ti{x),u=ui{x) 

For the non Abelian complexifier we find 

Zi(.) = ^i(x) - E 4±~L (2.29) 



while for the Abelian one we obtain 

Zi{x) = Aiix) - ^ E ^ E.ipi^ixMx)) (2-30) 

Notice that in both cases the imaginary part of Zi is only quasi local in Ej, that is, we can recover from 
Zi only up to the resolution provided by the parquettes. 
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2.3 Coherent States for Background Independent Gauge Theories 

We now come to compute the coherent states. The first step is to write the 5 distribution as 

5a, = Y, ^«(^o) <^^'- > (2-31) 

s 

where the sum is over ah spin network labels s = {'y,TT,m,n). Hence the coherent state is given by 

V^z = 5^ T,(Z) <e~^/^r„. > (2.32) 

s 

Here C is obtained by replacing in (|2.19|) . (|2.26|) or (|2.27|) respectively the gauge covariant flux, area or flux 
functionals by the gauge covariant flux, area or flux operator [IH [HI HI] respectively which are positive, self 
adjoint operators with pure point spectrum only. 

It remains to compute the action of the heat kernel and for this purpose we restrict to the case D = 3 
of ultimate interest. Again we do this separately for the two types of comp lexifiers. 

2.3.1 Gauge Covariant Flux Coherent States 

There is in principle an operator ordering problem involved in the quantization of (j2.19p . however, the 
regularization of [2j shows that there is no action of the operator valued distribution *E{x) on a holonomy 
A{p) if *E(x) is smeared over an infinitesimal surface element of a surface in which the path p lies. Let us 
introduce the matrices 

Ojkia) ■■= -\ T,{Tkkdg{Tj)) (2.33) 

where we have assumed the normalization Tr(rjrfe) = —25jk)- Since G is compact, we can always embed 
into a subgroup of some U{N) so that tJ = —tj, = whence Ojk{g) is real valued. Moreover, the 
obvious identity Ojk{g) = Okj{g^^) as well as the fact that Ad acts on Lie(G) whence Adg{Tj) = Ojk{g)Tf^ 
reveals that 

0jk{9)0ji{g) = 6ki (2.34) 

so that g i— > Ojk{g) is a subgroup of 0(dim(G)). 

The known quantization of the non gauge covariant flux [21 [JT] together with the above mentioned trivial 
action on Ojk{A{ps{x)) now reveal that 

E^)T^,j,m,n = i^l E '^x{S,e)0,k{A{ps{x)))hc^,T^,,,m,n (2.35) 

ee-E{7) x^Sr\e 

where is the right invariant vector field of G acting on g = A{e), specifically = Tt{Tjgd/dg'^). Here 
we have assumed that the graph has been adapted to S by suitable subdivisions of edges, such that each 
edge of 7 is either outgoing from an isolated intersection point or completely lies within S or lies completely 
outside S. 

Formula ()2.35p can now be plugged into (j2.20p . Since again there is no action of widehatE{S) on ps{x) 
we find 

^V^,i,m,n = -4 E E ^-(^'^) E ^y{S,e)Oui{A{ps{x)-^ops{y)))^X^,Xi,T,^,,m,n 

e,e'e-E(7) xeSne j/e5ne' 

(2.36) 

The appearance of the matrix Oki{A[ps{x)~^^ o ps{y))) makes the computation of the spectrum of (|2.36p 
rather difficult for a general graph. However, it becomes simple in case that the graph is such that the 
surface S has only a single isolated intersection point x with the graph. In that case (|2.36p becomes 

Ejisfr^j^^^n = -iU Yl E ^-('^' ^) i^e? T^,J,m,n (2.37) 
ee-B(7) xeSne 
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One can now introduce similar as in |41j the vector fields 

yj± = -i xi/^^n = yt+Yt (2-38) 

cr,(e,S)=±l 

so that we obtain the linear combinations of Casimir operators 

^) = f [2(yJ+)2 + 2{Yff - [Y^f] r,,,- (2.39) 

A special case arises when x is an interior point of a single edge e = (ei)~^ o 62 intersected transversely so 
that axiS,ei) = — cra;(5, 62) = ±1 and thus T^j,m,n is gauge invariant at x. Then ()2.39l further simplifies to 

Ej{sfT^j,m,n = 4 [-iXi/2f T^j,m,n (2.40) 

For G = U{1)'^ or G = SU{2) respectively the eigenvalues of (— iX|)^ are given by (ni)'^ and je{je + 1) 
respectively. This special case arises for the case of the polyhedronal cell complex complexifier when 7 
is a graph dual to it, that is, there is precisely one edge e of 7 which intersects a given face S and if so 
transversely. 



2.3.2 Area Coherent States 

Notice that for each direction I, each graph 7 and each stack a the Lebesgue measure of the set of times 
t such that contains a vertex of 7 or that p^j contains entire segments of edges of 7 vanishes. By the 
properties of the area operator and flux operator, it follows that those points do not contribute to the heat 
kernel evolution and therefore we may assume without loss of generality that each p^^ intersects the edges 
of 7 at most transversely in an interior point. Now consider in the non Abelian case for natural numbers 
A^e e No the set 

:= {t G M; \pi, H e\ = N^} (2.41) 

where we have abbreviated N := {Ae}egE(7)- This is the set of parquettes within stack which intersect 
edge e precisely times transversely. Likewise, consider for the Abelian case for integers Ng G Z 

5^°^:={teM; Yl <^{pit,e)p = N,} (2.42) 

where for any surface S intersecting e transversely, the number (t{S, e)p for p £ S De takes the value +1 or 
— 1 respectively if the orientations of S and e at p agree or disagree respectively. This is the set of parquettes 
within stack p^ whose signed intersection number with edge e is precisely N^. 
In both cases let 

:= / dt (2.43) 

■J 

be the Lebesgue measure or length of those sets. These length functions are needed in order to define a 
cylindrically consistent family of heat kernels as was first observed in [38] . Then the action of the complexifier 
on SNWF's is diagonal 

jTs = XsTs (2.44) 
The corresponding eigenvalues are given for G = SU (2) by 

= 1^ E E [ E Vie(je + 1)]^ (2.45) 

I,a N e£E(-y) 
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while for G = C^(l)'^ they are given by 

= ^ E E [ E nif (2.46) 

-f.a N eG-E{7),i 

Here we have used that the irreducible, non trivial representations of SU (2) are given by positive, half integral 
spin quantum numbers je 7^ while for U{Vf' they are given by triples of integers rvi 7^ 0, j = 1,2,3. 
Furthermore, with k = SvrGNcwtoni = is the Planck area. The ratio t := i'j,/a'^ is known as the 
classicality parameter. Without dynamical input, this is a free parameter for our coherent states that 
decides up to which scale the fluctuations of operators are negligible. 



2.4 Cut - Off Coherent States 

Formulae (I2.32p . (j2.36p ()2.45p and ()2.46p display the coherent states in closed form. Unfortunately, although 
the eigenvalues of the heat kernel grow quadratically with the representation weight, these states are still 
not normalizable because the Hilbert space is not separable, or in other words, the SNWF's are labelled by 
the continuous parameter 7. In view of the uniqueness result when insisting on background independence, 
the non separability is not avoidable and one must accept it. The observation is that ()2.32p defines a well 
defined distribution on the dense subset of 7i consisting of the finite linear span of SNWF's. To extract 
normalizable information from ipz we introduce the notion of a cut - off state labelled by a graph 7. These 
are defined by 

i;z,^:= ^ Ts{Z) < e-^/^T,, . > (2.47) 

That is, the sum over all spin networks s = (7(5), 7r(s), m(s), n(s)) is truncated or cut off to those whose 
graph entry 7(5) is a subgraph of the given 7. The Ansatz is then to use ipz,'y for suitable 7 as a semiclassical 
state. 

Notice that both (j2.45p and (j2.46p respectively can be rewritten in the form 

= ^ E C,e' VjeUe + l) VJe'Ue' + l) (2.48) 



2 

e,e' 



and 



where the edge metric 



IE Q,e'r^i<' (2-49) 



2 

e,e' 



C,e' = E E ^n' ^^^e' (2.50) 
I,a N 

has entered the stage. Such non diagonal edge metrics have already appeared in other background dependent 
contexts [Ml ES] . The edge metric decays quickly off the diagonal because for most edge pairs e / e' there is 
no direction and no stack in that direction intersecting both e, e' which means that = for N^, N^, / 
for such edge pairs. It is for this reason that we will be able to actually carry out our calculations. 

Using the edge metric, formulas (j2.29p . (j2.30p and (j2.45p . (j2.46p admit an interesting reformulation: 
The signed intersection number between a path e and a surfaces S is defined by (adopting convenient 
parametrization) 



a{S,e) := / dx" / dy" dy^ - eab^ 6{x,y) = dt d'u [eabce"(t) 6{e{t),S{u)) 

Je Js 2 Jo 7[o,i]2 du^ du^ 

= Yl ^-('^'^) (2.51) 



xesne 
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while the intersection number is given by 

\a\{S,e) :== dtj^^^^^ <fu \e,,^e {t)^^^^^^\ 6{e{t),S{u)) (2.52) 

Both expressions can be regularized in such a way that entire segments of e that lie inside S do not contribute 
to the integral [H]. Notice that |(T|(e,5) / |(T(e,5)|. Then it is not difficult to see that for SU{2) 

^le' = Y. I dt\a\{e,pf)\a\{e',pf) (2.53) 

while for U{1)^ 

^le' = Yl J dta{e,pf)a{e',pf) (2.54) 

a,I 

To verify (j2.53p . (j2.54p it is easiest to use directly the action of non Abelian area and Abelian flux operators 
respectively on the corresponding SNWF [41j (with only transverse intersections) 

AriS)T^j,m,n = kl(e,5) x/mTl)] 

ee£{7) 

Ej{S)T^,n = ^llYl ^(e>'5)ni]r^,„ (2.55) 

and to plug this formula into the expression for C. An alternative proof is by realizing that in the non 
Abelian or Abelian case respectively 

ee-E(7) ee-B(7) 

where xs denotes the characteristic function of a set. When plugging (|2.56p into (j2.50p and solving the 
Kronecker 6's when carrying out the sum over the integers N, (j2.53p and (j2.54p respectively result. 
From the easily verifiable properties of the (signed) intersection numbers 

cr(eoe', S) = a{e, S)+a{e' , S), cr{e~^ , S) = —cr{e, S); \a\{eoe' , S) = \o-\{e, S)+\a\{e' , S), \a\{e~^ , S) = \(j\{e, S) 

(2.57) 

it follows immediately that 

P{e oe',eo e') = P{e, e) + P{e', e') + 2V{e, e'), r(e-\ e~^) = V{e, e) (2.58) 

This is precisely the generalization to non diagonal edge metrics of the cylindrical consistency conditions of 
the complexifier |22i [38] . Notice that for the general area complexifier (j2.24p we arrive instead at the edge 
metrics 

,= / dii{S) I d^l{S')\a\{S,e) K{S,S')\a\{S' ,e'), ll^, = f dn{S) f d^iS) a{S,e) K{S, S') a{S' ,e') 

(2.59) 

Finally we have for any edge e 

f f /■! (n^ E^)(e(s)) r 

/ dx'^ta'[Zi-A^^]{x)=Y, dtAripi-) ds / d'u[e%s)nf\u)6{pf{u\e{s))] 
Je t^J Jo ^[{niE%e(sW J 

(2.60) 

in the non Abelian case while for the Abelian case 

/ dx" ia^Zi - A^^]{x) = Y. [ dt Ejipf) a{pf,e) (2.61) 

Interestingly, if E does not vary too much on the scale of a plaquette, then ()2.60p actually reduces to (I2.6ip 
which is written directly in terms of the signed intersection number and plaquette fluxes. This will be useful 
later on when we compute expectation values. 



/ 
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2.5 Comparison of Gauge Covariant Flux and Area Coherent States 

Consider the case that the plaquttes are much smaller than the edges with respect to the three metric to be 
approximated by the coherent states and that the edges do not wiggle much on the scale of the plaquettes. 
Then for each I the number of stacks that do not contain a vertex of 7 but still intersect 7 drastically 
outnumbers the number of stacks that do contain a vertex. Moreover, among the edge free stacks, the 
number of stacks that intersect only one edge completely outnumbers the ones that intersect more than one 
edge. Finally, among those with single edge intersections, the number of stacks that intersect the respective 
edge once completely outnumbers the ones that do more than once. For this reason, in these cases the 
expressions (I2.45P and p.46p can be replaced with good approximation by simpler expressions of the form 

^^ = ^ E l2jeUe + l) (2.62) 

and 

^^ = ^ E ^2KF (2.63) 

eeB(7) 

respectively where the length function Q = g2e solves Ij^^, = le + le', K-^ = in order that the complexifier 
has cylindrically consistent projections. This is the form of the heat kernel eigenvalue considered for the 
states in |38) . As shown in |22] , these eigenvalues cannot come from a known classical complexifier so that the 
complexification map A ^ Z, without which the Z label of the coherent state has no relation to the phase 
space point to be approximated, is unknown. When using the complexifier coming from a polyhedronal cell 
complex, a concrete relation between Z and the phase space can be given for specific graphs, the above 
eigenvalues arise as we saw in section 12.3.11 and as shown in [22j . 

Let us also check that the area complexification map Z in (j2.29p and ()2.30p comes close to the gauge 
covariant flux one ()2.22p . at least on certain graphs. Let 7 be a graph dual to the cell complex P. Thus, 
for each edge e there is a unique face Se which intersects e in an interior point transversely such that 
(Js^ne{Se, e) = +1 and no other face intersects e. Then the gauge covariant flux complexification map at the 
level of the holonomies is given by [llj 

A{e) ^ ge{Z) := Z^(e) = exp(-irj^j(5e)/L2) A{e) (2.64) 

For SU{2), iTj are the Pauli matrices while for ^7(1)'^ iTj = 1, j = 1,2,3. In contrast, the area complexifi- 
cation map is given at the level of the holonomies by 

A{e) ^ Z{e) = Pexp( / Z^tj) (2.65) 



where V denotes path ordering where Zi is given in (12.29P and (|2.30p respectively. Now for sufficiently 
"short" edges we have Z{e) ~ exp(J^[Z — AyTj)A{e) to leading order in the edge parameter length. If we 



assume that Ej is slowly varying at the scale of the plaquettes then we have ^'^iPaj(^x)ti{x)^ ~ y [-^j (^)"'a(^)]^ 
so that (|2.29p is approximated by 



{Z^ - A^) - E / ^^7T7^ nl{e{t))E^M^)) (2.66) 

where we have assumed that e is the embedded interval [0, 1]. Now consider the case that the graph is in 
fact cubic and that the stack family and the graph are aligned in the following sense: 

Suppose that we have an embedding X : — > u; s 1— > X{s). For ejjx = 1 we define Xl{u^,u^) := 
X{s^ = t, s'^ = u^,s^ = v?). This defines linearly independent foliations with leaves Ljt = Xj-i^). 
The corresponding stack families are labelled by a := (a^,a^) S 1? and defined by : [0,1)^ — > 



17 



cr; X^^{u) := Xl{[a^ + u-*^]/, [a^ + u'^]l) where / > is a certain parameter. The edges of the cubic 
graph are labelled by verticies v = G Z'^ and directions / and are defined for ejjx = 1 by 

Cyj : [0, 1] — > a; e^j{t) := X(s^ = [v^ + 1]6, s"^ = v^6, = v^5) where (5 > is another parameter. 
In this situation, (j2.66p can be further simplified to 

r (z^- _ A^) « _ r ni{e,j{t))E';iej,it)) ^ -^5E^{pi) (2.67) 
Je^l 0. Jo a 

where pi is any plaquette in the stack in / direction intersected by ejy. 

Thus, for cubic graphs, which are the only ones considered so far in semiclassical calculations, we get a 
close match between (|2.64p and (|2.67p whenever the cubic graph and the stack families are aligned. However, 
there is still an important difference: 

The parameter area P of the plaquette pi in (j2.67p has no a priori relation to the parameter length 6 of the 
edge Syi while the parameter area of the dual face Se.^j in (j2.64p is of the order 5^. These considerations 
reveal that the individual plaquettes of the stacks cannot be interpreted as the faces of a dual graph although 
roughly [5/1]'^ of them combine to a face. Hence the states considered in [22| are genuinely different from 
those in [TT] . 

This will turn out to be important: 
We will see that in order to be able to perform practical calculations for SU{2) with off diagonal edge 
metrics, we need i <C 5 in order that the edge metric is close to diagonal for generic graphs. It turns out 
that if we use the same parameter a in the label Z of the state and for the classicality parameter t = I'^p/a^ 
then the expectation value of the volume turns out to be of the order of (1/6)^ too small. Hence, there is a 
tension between the possibility to perform practical calculations and the correctness of the classical limit. 
The only analytical calculation possible with / = 6 uses a graph which is aligned with the stacks and thus is 
necessarily cubic. While the result of that calculation results in the correct classical limit, this calculation 
is of limited interest because we saw already above that for this case the coherent states of [22] reduce to 
those of [11] for which we knew already that the classical limit is correct. 

However, the purpose of this paper is to test the semiclassical limit for graphs of non cubic topology. 
This can be done with the states of [11] without limitation. With the states of |22j this is possible if we 
redefine Zi ^ Aa + frC-^a — Aa) where b <^ a. This rescaling is actually not in the spirit of the complexifier 
programme, but it repairs the semiclassical limit of all operators built from the fluxes. It will then turn out 
that for graphs that satisfy 1/6 = b/a the correct classical limit results for n = 6 only. As already mentioned 
in the introduction, one could rescale the label of the coherent state by a different amount in order to reach 
the correct semiclassical limit of the volume operator for one and only one n ^ 6. However, that would 
destroy the correct semiclassical limit of other operators such as areas. Hence the rescaling by is 
harmless in the sense that it reproduces the semiclassical limit of all operators while n— dependent rescaling 
do not. 

Also with respect to the states of [11] the value n = 6 is singled out. The fact that the cut off states of 
|22j have acceptable semiclassical behavior only when the corresponding cut off graph and the label of the 
coherent state satisfy certain restrictions imposed by the structure that defines the complexifier, in this case 
the size of the parquettes, is similar to the restrictions imposed on the by the polyhedronal cell complex 
complexifier [llj, namely that the graph be dual to it. 

2.6 Justification for Replacing SU{2) by U{1)^ 

The considerations above have revealed that practically useful cut - off states will be based on graphs which 
are much coarser than the parquets so that the edge metric is diagonal in very good approximation. We will 
restrict to such graphs in the calculations that follow and find independent confirmation for that restriction 
as well in the form of the quality of the semiclassical approximation. Assuming exact diagonality and thus 
suppressing the corrections coming from off - diagonality which we will show to be small under the made 
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coarseness assumptions, the cut - off states in fact factorizes 



where for 5C/(2) 



while for C/(l)" 



'^Zn= n ^^'7,e (2.68) 
ee-E(7) 



V^z,,,e(^) = + 1) e-t'^^(^+i) X,{9eA{e)-^) (2.69) 

2j=0 

V'z,7,e(A) = e-5'^E,(n^)^ Xn{9eA{e)-^) (2.70) 



Here Xj and Xn respectively denotes the character of the j— th and n— th irreducible representation of SU{2) 
and U{1)^ respectively. 

Under the made assumptions, the edge metrics Q are identical for both groups because, while is 
defined for non negative integers N only in the case of SU{2) while for C/(l)^ all integers are allowed, for 
the graphs under consideration for each edge e only either N(, = +1 or iVg = — 1 leads to non vanishing l^j^^ 
so that these numbers in fact coincide and since we take the diagonal elements of the edge metric (|2.50p 
both signs lead to the same Q. 

Finally, if {Aq,Eq) is the phase space point to be approximated and from which we calculate Z = 
Z{Ao,Ao) via ([2:^9]) and (12:3(1 then for SU{2) we have, 

<7e « exp(-ir,/^(e)) e^Tj Ja,), P^ie) = ^ ^ '^^^mT K(<t))ni{em (2.71) 

while for U{1)^ we have 

9e = {gi)U^ 9i = exp(-/^(e) + i f J^^ (2.72) 

J e 

where as before we have made the approximation 



which is valid if Eq is slowly varying at the scale of the plaquettes. 

Thus, given Z = Z{Aq,Eq), we have the following abstract situation under the made assumptions: 
1. 

For each edge e there exist vectors Pq (e), ^Q(e) such that for SU (2) we have ge ~ exp(— irjPjrj) exp(rjAQ(e)) € 
SL{2,C) = SU{2f while for U{lf we have ge = (e-^(^)+^^o(e))3^^ G (C - {O})^ = {U{lff. 2. 
The coherent states adopt approximately product form il^z,'y ~ nee-E;(7) V'g^ where 

00 

V',(/i) = + 1) e-«^^(^+i)/2 xM-') (2.74) 

2j=0 

for h G SU{2) while 

^g{h) = Y e-*'^^?-"? Xnigh-') (2.75) 

for h G C/(l)3. 
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Now, as anticipated in the introduction, using the tools of semiclassical perturbation theory [ITO] we are 
able to calculate to expectation value of the volume operator V of LQG with respect to the correct SU{2) 
coherent states in terms of the expectation value of of a certain operator Q, where V = which we 

display explicitly in the next section and which is a sixth order polynomial in the right invariant vector 
fields Xe on SU{2) where Xi acts on /ig in (j2.74p . The crucial observation, made in [11], is that if we 
simply replace the SU{2) right invariant vector fields in Q by U{1)^ right invariant vector fields Xl acting 
on /ig in (|2.75p and if we replace the SU{2) coherent states ()2.74p by the related U{1)'^ coherent states in 
(j2.75p . then the remarkable fact is that the expectation values of polynomials of right invariant vector fields 
actually coincide to zeroth order in h. By the same argument, this will be also true if we perform the right 
invariant vector field replacement already at the level of V rather than Q. This observation was also key in 
the semiclassical analysis of \20\ [21] . 

This feature is maybe not as surprising as it looks at first sight because, after all, the coherent states 
for both groups are to approximate the same phase space points. The underlying reason is that that the 
classical phase space of the SU{2) theory (i.e. the range of fields and the symplectic structure) and of the 
fictive U{lY theory actually coincide. It is only when we add the dynamics of the theory as for instance 
the Gauss constraint that we see a difference. The Gauss law is taken into account in two ways, first by 
using the appropriate group coherent states, here SU{2) or f7(l)^ respectively, which is dictated by the fact 
that the underlying holonomies take values in the appropriate group. Secondly, one can construct quantum 
Gauss constraint invariant coherent states [111 B2] by averaging over the gauge group action at the verticies. 
Denote this group averaging map by r]. Then, as shown in [111 142] . we have that < rj{'ipz,'y, ^'ni'^z,'y) > and 
iljz,'y, ^■0z,7 > agree to zeroth order in h (notice that the Gauss invariant Hilbert space is an honest subspace 
of the kinematical Hilbert space so that the same inner product can be used) for any Gauss invariant opera- 
tor A such as the volume operator because the overlap function between coherent states peaked at different 
phase space points is sharply peakec|l|. This justifies the use of the kinematical states employed in this paper. 

To summaries: 

Using kinematical C/(l)^ coherent states is a convenient approximation for actual SU{2) coherent state ex- 
pectation value calculations for Gauss invariant operators if one is only interested in the zeroth order in h. 
At non vanishing orders in h there will be differences but we are not interested in them for the purposes 
of this paper. One may wonder whether the argument made above, namely using kinematical rather than 
Gauss invariant coherent states also survives when considering the spatial diffeomorphism constraint. This 
issue, currently under investigation, is more complicated in part because it is not completely obvious which 
distributional extension of the classical diffeomorphism group one should use [l3]. For the purposes of this 
paper this is of no concern because we are looking at the local volume operator which is not spatially diffeo- 
morphism invariant so that expectation value calculations with respect to spatially diffeomorphism invariant 
coherent states are meaningless. It is the local volume which enters the Hamiltonian and Master constraint 
and verifying the semiclassical limit of those only makes sense at the kinematical Hilbert space level (one 

^In more detail we have 



ri{ipzn) = 




(2.76) 



where ag{^z.-^){A) = ipz,-y{ctg{A)) and [Qg(A)](e) — g{b{e))A{e)g{f{e)) ^ where b(e) and /(e) respectively denote beginning 
and final point of e respectively. Now due to gauge covariance of the coherent states we have ag{tpz,-y) ~ tpa _i {z),-y so that the 
gauge invariant coherent state expectation value of a gauge invariant operator becomes (using the invariance properties of the 
Haar measure) 

< v{'^z,i),Ari{ipz,^) > ^ JglvMl Ylvevji) ^Mg(gt.) < i^a,(z),i, Aipz,^ > 

\\v{i'Z,-YW /gIV(t)I Uv€V(-y) '^^^H(9v) < ^ag{Z),7.'^Z,7 > 

Now from [11] we know for gauge invariant polynomials A in right invariant vector fields that the peakedness property 

< i'z',„Ai,z,-, ^'^I'lT'^tf^ ^ < i-z'.^^z., > [1 + Om (2.78) 

\\YZ,-y\\ 

holds. Now the claim is immediate. 
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cannot check the correct classical limit of a constraint on its kernel). Once this limit is verified, one has 
confidence that the physical Hilbert space defined by the Hamiltonian constraint is correct. 

3 Regular Simplicial, Cubical and Octahedronal Cell Complexes 

In order to perform the calculations in our companion paper for the coherent states of [52] , we need the 
specific embedding of the n = 4, 6, 8 valent graph relative to the stack families. This can easiest be done 
by starting from regular dual simplicial (tetrahedronal) , cubical and octahedronal partitions of the three 
manifold a. For the coherent states of \lT\ this section is not needed except that it shows the existence of 
(regular) polyhedral cell complexes dual to n = 4, 6, 8 valent graphs such that all cells of that complex are 
platonic solid bodies, i.e. tetrahedral, cubes and octahedra respectively. 

In fact, it is possible to define such partitions all from refinements of cubical decompositions such as 
sketched in figure [TJ We perform the analysis for each chart X : — > cr separately and use the Euclidean 




Figure 1: Cubic cell decomposition. 



metric on M.^ in the following definitions. 

Definition 3.1. 
i. 

A cubical partition ofM.^ is defined by the cubes Cn, n^H? where 

Cn = {s £ / = + t^, 1=1,2, 3} (3.1) 

The boundary faces (squares) of Cn are taken with outward orientation, 
ii. 

A simplicial partition of M.^ subordinate to a cubical one is defined as follows: 

First draw in C(o^o,o) diagonals on the boundary squares such that the diagonals on opposite squares are 
orthogonal. Specifically, in the face defined by s^ = 0; s'^,s^ G ^ijk = 1 th^ diagonal is the line 

t {s^ = 0,s-^ = t,s^ = t), t G [0, 1] while in the face defined by s^ = 1; s'^,s^ £ [0, 1]^; cijk = 1 the 
diagonal is the line t ^ {s^ = 1, s"^ = t, = 1 — t), t G [0, 1] ■ 

Now continue this pattern of orthogonal diagonals in opposite faces to the six cubes adjacent to cq where 
common faces have the same diagonal. This also defines the remaining four diagonals in those six cubes by 
connecting the endpoints of the already present two diagonals. 
Finally continue this process for all cubes. 

The face diagonals define altogether five tetrahedra that partition each cube. We will take their boundary 

triangles with outgoing orientation. 

Hi. 
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An octahedronal partition of subordinate to a cubical one is defined as follows: 

For each cube draw the unique four space diagonals. Specifically in 0(^0 0^-^ these are the lines t {t,t,t), {t,t,l — 
t), {t, 1 — t,t), (1 — t,t,t); t £ [0, 1]. These partition each cube into six pyramids with common tip in the 
barycenter of the cube and with the six faces of the cube as their bases. Now glue two pyramids in adjacent 
cubes along their common base. Obviously, two glued pyramids define an octahedron which we take with 
outgoing orientation. 

The basic building blocks of the tetrahedronal and octahedronal decompositions are displayed in figures 
[21 [3] and U] respectively. When gluing the bases of the pyramids along the faces of the original cubes one 




Figure 2: Type A triangulation of a cube. 




Figure 3: Type B triangulation of a cube. 

obtains an octahedronal decomposition as displayed in figure O It is maybe not completely obvious that the 
drawing of the diagonals that are to define the tetrahedra is a consistent and unique prescription. To see this, 
we use the checkerboard visualization displayed in figurefig7: First draw all plaquettes in the s"^ = n £ Z, 
layers. Now take the n = layer and draw the diagonal for the plaquette in that layer that belongs to 
C(o,o,o) &s prescribed in the definition. Define that plaquette as "black". Now turn the n = layer into a 
checkerboard in the unique way consisting of black and white plaquettes. The other layers n 7^ are also 
turned uniquely into checkerboards by asking that checkerboards in adjacent layers are complementary, i.e. if 
the plaquette {n^,n'^,n^) is white (black) then the plaquette (n^, n^, it 1) is black (white). Draw diagonals 
in plaquettes of opposite colour orthogonally to each other. This defines face diagonals in the s^ = nconst. 
layers. These have the property that they form squares in each layer which lie at an angle of 7r/4 relative to 
the plaquettes and which are such that only every second plaquette corner is a vertex of these squares. We 
will refer to such corners that are verticies as "used". It is easy to see that in adjacent layers, used plaquette 
corners lie above unused ones. Now draw the remaining face diagonals in the s^,s'^ = n =const. layers by 
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connecting the used corners in adjacent layers using the appropriate diagonals of the cubes. This results 
in the triangulation depicted in figure [71 We now define the graphs dual to these particular polyhedronal 
decompositions. 

Definition 3.2. 

The graph in dual to the above simplicial, cubical and octahedronal cell complexes is obtained by connecting 
the barycentres of adjacent tetrahedra, cubes and octahedra respectively by straight lines through their common 
triangles, squares and triangles respectively. Here the barycenter of a region R CB. is defined as usual by 

B{R) = (3.2) 
•J R 

The advantage of the explicit definition of the cell complex is that we can explicitly label the edges and 
verticies of the dual graph. This is of course only feasible for sufficiently regular graphs, otherwise we run 
into difficult bookkeeping problems. 

1. Cubical Graph 

The barycentres of the cubes c„ are evidently the points w„ := {n^ + |,n^ + ^,71^ + 1) which form the 
verticies of the dual graph. The edg es 6ri7) I — 1)2,3 which connect the verticies with labels n and 
n + bj respectively, where bj is the standard unit vector (bj)'^ = 5j, have the explicit parametrization 
e-n,i{t) = Vn + tbi, t S [0, 1]. The other three edges adjacent to Vn are ingoing and are given by Cn-bij- 
These edges form the 1 skeleton of another cubical cell complex which is just shifted by the vector 
(|, |, ^) from the original one. 
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Figure 6: Checkerboard visuahzation of the triangulation. 




Figure 7: Triangulation. 



2. Tetrahedronal graph 

The tetrahedronal graph is the most complicated one because there are two different types of simplicial 
decompositions of a cube into five tetrahedra. Type A. corresponds to the case that the verticies of 
the internal tetrahedron within a standard unit cube are given by (0, 0, 0), (1, 1, 0), (1, 0, 1), (0, 1, 1) 
while type B. has verticies at (1,0,0), (0,1,0), (0,0,1), (1,1,1). These types alternate in adjacent 
cubes as we move in any of the three coordinate directions. Hence, by defining the cube cq to be of 
type A., the triangulation is completely completely specified. Indeed, the type of Cn is A. if + n? + 
is even and of type B. otherwise. 

To determine the dual graph, we first discuss the barycentres of the tetrahedra for the two types 
separately for a standard unit cube as well as the edges of the dual graph that lie within it. The 
verticies of and the edg es in Cn follow then by translation by n — n bj. Notice that a tetrahedron 
T based at v and spanned by vectors ej, that is, T = {v + t^ej; < < 1; + t"^ + < 1}, has 
barycenter at B{T) = v + |(ei + 62 + 83). 

A. Type A. 

The barycenter of the interior tetrahedron coincides with the barycenter vq := |(1,1,1) of the 
cube. The barycentres of the remaining four exterior tetrahedra based at verticies 
(1,0,0), (0,1,0), (0,0,1), (1,1,1) respectively are at < := i(3,l,l), = i(l,3,l), = 
j(l,l,3), = ;j(3, 3,3) respectively. Accordingly, the dual edges within the cube are = 
- Vo, a = 1,2,3,4. 

B. Type B. 

The barycenter of the interior tetrahedron coincides with the barycenter vq := ^(1, 1, 1) of the 



24 



cube. The barycentres of the remaining four exterior tetrahedra based at verticies (0, 0, 0), (1,1,0), (1,0,1), (C 
respectively are at := |(1,1,1), = |(3,3,1), = ;j(3, 1,3), vf = |(1,3,3) respectively. 
Accordingly, the dual edges within the cube are = — vq, a = 1, 2, 3, 4. 

It remains to describe the dual edges that result from gluing the faces of the exterior tetrahedra of 
adjacent cubes. But this is simple because each of the exterior tetrahedra within a cube has three 
triangles as faces which lie in the three coordinate planes, hence the gluing is between those triangles 
which result from drawing the respective face diagonal within a boundary square of a cube. Hence, 
each cube has twelve edges perpendicular to the twelve boundary triangles of the exterior tetrahedra 
which are adjacent to the four barycentres of those exterior tetrahedra. Altogether we can identify 
six possible gluings, namely either going from type A. to type B. when moving along the positive / 
direction and gluing along the =const. plane or going from type B. to type A. when moving along 
the positive / direction and gluing along the =const. plane. As one may check, the type A. to type 
B. gluing in / direction corresponds to two dual edges running from vf to a and from to vf 
types of verticies respectively. Likewise, the type B. to type A. gluing in / direction corresponds to 
two dual edges running from Vj to and from to Vj types of verticies respectively where ejjK- 
In all cases, these I direction edges have coordinate length ^ as one may easily calculate. 

Altogether, we can now easily describe the dual lattice as follows: 

The verticies are labelled Vn,ai « = 0, 1, 2, 3, 4 with Vnfi = n -\- vq and Vn,a = n + v^, a = 1,2,3,4 if 
v} -\- v? -\- is even while Vn,a 

= n + B^, a = 1, 2, 3, 4 if + + is odd. The edges are labelled by 
en,a, a = 1, 2, 3, 4 and en,/j, 1 = 1,2,3, j = 1,2 where en,a{t) = n+vo+t{v^— vq if n^+n'^+n^ is even, 
<in,a{t) = n + VQ + t{v^ - Vq) if n^+n^ + n^ is odd, en,i,i{t) = n + vf + ^hi and enj,2{t) = n + vf + ^bj 
if + 71^ + is even and finally en,i,i{t) = n + vf + ^bi and en,i,2{t) = n + vj^ + if + + 
is odd where eijK = 1- 

3. Octahedronal Graph 

Each cube contains six pyramids or halves of the octahedra. Therefore the barycenter of an octahedron 
coincides with the barycenter of the common boundary face of the two cubes that contain it. It follows 
that the octahedra may be labelled by Onj corresponding to the verticies Vnj = n+^bj+^bx', ^ijk = 1 
which define its barycenter. Such an octahedron has the property that it has a common base of two 
pyramid halves which lies in the =const. plane. For the vertex Vnj we define four edges en,i,.j,(7, J / 
/; (7 = lb outgoing from it through the explicit parametrization enj,j{t) := Vn,i + ^{bi + crbj) which 
connects the verticies Vn,i and v^^i^-^^j^^^i^^ j. Notice that these edges lie in the {I, J) or {I,K) plane 
but there are no edges in the (J, K) plane adjacent to Vn,i- The other four edges adjacent to Vn,i have 
ingoing orientation. 

As an aside, notice that the 1-skeleton of an octahedral cell complex as defined above is an eight valent 
graph after removing the edges of the original cubes. 

The basic building blocks of the dual graphs are displayed in figures [8l O [10] and [11] respectively. The 
connection of the tetrahedronal lattice with the diamond lattice is as follows: 

For each cube of type A. or B. respectively, keep the interior tetrahedron. Now move the barycentres of the 
remaining exterior tetrahedra into that corner of the cube which is also a corner of the tetrahedron under 
consideration. In this process, the edges dual to the faces of the interior tetrahedron become halves of the 
spatial diagonals of the cube. Finally drop all the other edges which were running between the barycentres 
of the exterior tetrahedra. The result is a diamond lattice. Its basic building blocks are depicted in figures 
[T2] and [13] respectively. It is also four valent, however, it does not have a piecewise linear polyhedronal 
complex dual to it (i.e. whose faces (which are subsets of linear planes) are in one to one correspondence 
with the edges) . It does have a cell complex dual to it if one gives up piecewise linearity by suitably rounding 
off corners but that is inconvenient to describe analytically. On the other hand, the natural polyhedronal 
complex consisting of the interior tetrahedra of the original cubes with the cubes deleted consists of those 
tetrahedra as well octahedra which surround half of the corners of the original cubes. Only half of the 
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triangle faces of those octahedra are penetrated by the edges of the diamond lattice. The building of this 
semi dual polyhedronal cell complex consisting of tetrahedra and octahedra respectively is visualized in 
figures [HI [I5l [El [17] and [18] respectively. In order to achieve the desired duality, one has to fill in the 
original cubes again which then triangulate those octahedra into eight tetrahedra. This then results in the 
additional verticies and edges that we described and depicted in figures [10] and [TTJ 

4 Volume Operator Expectation Values for Dual Cell Complex Coherent 
States 

In this section we compute the expectation value of the volume operator with respect to the dual cell complex 
coherent states of jll] . This section is subdivided into two parts. In the first we review the definition of 
the Volume Operator. In the second, we perform the actual calculation. In order to carry it out explicitly, 
we have to specify the graph and the dual cell complex. Here we focus our attention on arbitrary graphs 
with the following properties: 1. All verticies have constant valence n = 4, 6, 8 and 2. the dual cell complex 
consists only of tetrahedra, cubes and octahedra respectively. Such graphs and dual cell complexes exist as 
we showed explicitly in section [3] This is all we need for the purposes of this section, more specifics about 
the graph and the complex are not needed. 
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Figure 10: Type A triangulation of a cube and dual four valent graph. 



/ 



Figure 11: Type B triangulation of a cube and dual four valent graph. 



4.1 Review of the Volume Operator 

The classical expression for the volume of a region i? of a semianalytical three dimensional manifold a is: 



Vr := / d^x v^det(g) = / ^/\^^lE\ 

JR Jr 



(4.1) 



where Qab is the three metric. The version of the volume operator [9] consistent with the triad quantization 
[lUj that enters the quantum dynamics [,4j has cylindrically consistent projections Vr_^ given by 



v&V{'y)nR 



where 



\ o — , 

\ ei,ej,eK,I<J<K<N\v&eir]ejr]eK 



(4.2) 



(4.3) 



Here N denotes the valence of the vertex, £p = hn is the Planck area, 



ei{v) 



Ti{[Tihi]^d/dhi) are 



right invariant vector on SU{2) acting on the holonomy hi := A{ei) (irj = aj are the Pauli matrices) and 
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Figure 12: Type A diamond cell with occupied lower, left front cube. 




Figure 13: Type B diamond cell with unoccupied lower, left front cube. 



e(e/, ej, Ci^-) is called the orientation function which is defined as follows: 

{1, iff e/, ej, e^are linearly independent at v and positively oriented 
— l,iff e/,ej,ei^are linearly independent at v and negatively oriented (4.4) 
0, iff e/, ej, e^are linearly dependent at v 

Here we take the convention that the edges at v have been taken with outgoing orientation, hence if in 7 the 
orientation of an edge e adjacent to v is actually ingoing, just apply the above expression to ip' {..,h~^ , ...) := 

V'(..,/le, ..). 

From (j4.2p . we deduce that the volume operator is a sum of contributions, one for each vertex. Therefore, 
in the expectation value calculations that follow it will be sufficient to calculate the expectation values for 
each V^^v separately and then to add the contributions. Notice that each of these contributions is of the 
form V^^v = \/Q'y,v where Q^^v is minus the square of the expression appearing between the modulus labels 
|..| in (j3.2p and therefore is a sixth order polynomial in the SU{2) right invariant vector fields. 

We will now proceed to calculate the general expression for the expectation value of the volume operator 
for an n = 4, 6, 8 valent graph. 
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Figure 14: Dual diamond cell of type A. 




Figure 15: Dual diamond cell of type B. 



4.2 Expectation Values 

We can actually perform a full SU{2) calculation as follows: 
The coherent states are explicitly given by [llj 

oo 

^za= n ^z,e, ^zAA) = Y. ^~''^'^^^'^^i^9e{Z)A{e)-^) (4.5) 

ee£;{7) 2i=0 

where t = l\,/L?' and ge{Z) is given by (j2.64p . The volume operator expectation value is given by 

< V{R) >za= Yl < '^7,- >^,7 (4-6) 

v(^V{'y)r\R 

Notice that due to the product form (j4.5p . the expectation value < V-y^y >z,'y only involves the edges adjacent 
to V. Now as we saw in the previous section we have V-y^v = yQ-^- By the arguments presented in the 
introduction, the zeroth order in ^ of < V^^^ >z,7 is given by -y/< Q^^v >z,7- Since Q^^v is a polynomial in 
right invariant vector fields, the results of [11] reveal that that to zeroth order in h the expectation value of 
any polynomial in the right invariant vector fields ii'jjXi is simply obtained by replacing it by Ej{Se) which 
is given in (j2.19p . 
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Figure 16: Dual diamond cell of type B with only the central four valent vertex left. 




Figure 17: Dual diamond cell of type B with only the central four valent vertex left and keeping only the 
faces adjacent to the vertex. 



It follows that to zeroth order in h we have < Q^^v >z.7= [P'y,v{E)]'^ where 

P^^,{E) = ^ ee,e',e") e'''' Ej{Se) EkiSe') EiiSe") (4.7) 

ene'ne" 

Notice that for sufficiently fine graphs, we can drop the holonomies along the paths Pei-c) involved in the 
definition of Ej{Se) as we approach the continuum. It is then clear that the correct expectation value of the 
volume operator is reached provided that (|4.7|) approximates the volume, as specified by Ej, of the cell of 
the polyhedronal complex which is bounded by the faces Se involved in (j4.7p . 

To do this, we use the fact that for sufficiently fine graphs a polyhedron P in a dual to a vertex of the 
graph lies in the domain of a chart Y so that P is the image under y of a standard polyhedron Pq in M^. 
Introducing 

na{s) - 2 ^abc e (4.8) 
we immediately find with P = Y{Pq) 

Vol(P) = / cPx V|det(E)(x)| = / cPs \l\dei{E{s)\ (4.9) 

Jp JPo 

where 

EUs) = EUYis)) niis) (4.10) 
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Figure 18: Dual diamond cell of type B with only the central four valent vertex left and keeping only the 
faces adjacent to the vertex, highlighting its octahedronal cell structure. 



Now for sufficiently fine graphs (|4.10p is approximately constant over Pq so that 



where 



Vol(P) ^ \/\det{E{s)\^^^^^^ Volo(Po) (4.11) 



Volo(Po) = / dh (4.12) 



is the volume of the standard polyhedron with respect to the Euclidean metric on M^. 

The idea behind this rewriting is that the fluxes Ej{Se) can be approximated by specific linear combina- 
tions of the [Ej {s)]y (s)=v so that a direct comparison between (14. 7|) and ()4.12p is possible. This is because 
a boundary face S is also the image under y of a standard face in so that (dropping the holonomies 
along the pdx) as explained) 

Ej{S) = ^ ^ eabc dx' A dx' E-^ix) = ^UK ds^ A ds^ £:/(s) « [E] {s)]y is)=v ^'{S"^) (4.13) 

where 

F/(5°)= / ItuKds-^ ^ds'' (4.14) 
is the / component of the Euclidean flux through Sq. Thus, plugging (j4.14p into (14. 7p we find 



\P,AE)\"^ ^ V I det(^(5))^(^)^^ Volo(i;) (4.15) 



where 



Volo(t;)= /l^ E ^e,e',e")e^'^ Fj{S^,)Fj{SI,)Fk{SI„)\ (4.16) 

V ene'ne" 

It remains to compare ()4.12p and (|4.16p . All of this still holds for general graphs. We now specify to 
purely n = 4, 6, 8 valent graphs with the above specified properties in order to test the correctness of the 
expectation value for specific, simple situations. Thus we know that for each vertex v the faces Se dual to 
the edges e adjacent to v form the surface a tetrahedron, cube and octahedron respectively. Thus we just 
have to compare (14. 7p with the volume of such platonic bodies as measured by E^. We will discuss the three 
cases separately. 
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4.2.1 Tetrahedron 

A standard tetrahedron is the subset 

To = {s G : < < 1; / = 1, 2, 3, s'^ + + < 1} (4.17) 

It has four boundary triangles given by 

= {s G : / = 0, < s-^,s^ < 1, s"^ + < 1; e/jx = 1} 

tl = {s G : < s-^ < 1; / = 1,2,3, + 5^ + = 1} (4.18) 

We easily compute 

Volo(To) = I (4.19) 
o 

while (remembering that the surfaces carry outward orientation if the edges are outgoing from v) 

Fj{t'j) = \5jj, Fj{tl) = -\ (4.20) 
Let us label the edges adjacent to f by ei, .., 64 where is dual to Y{t^), j = 1, 2, 3, 4. Then 



Volo(^) = =J\l E ee,,e„e,) e^'^ Frif^^) Fj{tl) 

V l<j<k<l<4 

= ^ V|e(ei,e2,e3) - €(61,62,64) - £(61,63,64) - 6(62,63,63)! (4.21) 



which still depends on the sign factors. Hence the expectation value takes values in the range 0, |, 3 
none of which coincides with |. For the explicit four valent graph that we constructed in section ([3]) each 
triple among the four edges has linearly independent tangents at v and the expectation value is given by 

/2 1 

^ > g which is too large. 
4.2.2 Cube 

A standard cube is the subset 

Co = {s G : < < 1; / = 1,2,3} (4.22) 
It has six boundary squares given by 

s?+ = {s G : = 1, < s\ < 1; ejjx = 1} 

s?_ = {s G : = 0, < s\ < 1; ejjx = 1} (4.23) 

We easily compute 

Volo(To) = 1 (4.24) 
while (remembering that the surfaces carry outward orientation if the edges are outgoing from v) 

Fjis%) = a5ij (4.25) 

with (T = lb. 

Let us label the edge dual to Y{s^^) by eja- Then the expectation value becomes 



Vol(t;) = /|^ E e{eiai,eja2,eKa3) cric^2(y3 eijxl (4.26) 



'48 

I ,J,K;ai,a2,cr3 
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which again depends on the precise embedding of the graph. For an actual cubical graph constructed in 
section [3l the edges e/+ , e/_ are analytic continuations of each other so that the orientation factor vanishes 
if two or more edges carry the same direction label I. Otherwise there are more contributions. Which 
orientation factors are allowed has been analyzed in detail in [18]. In the case of the actual cubical graph 
we have £(6/^1, ejo-2, ei^o-s) = 0'icr2cr3e/j/^ so that (j4.26p becomes 



Vol(-) = jl^ E ^wl = l (4.27) 



'6 

I,J,K 



which coincides with (I4.24p . 

4.2.3 Octahedron 

A standard octahedron is the subset 



Oo = {5 G M3 : \s^\ < 1 \s\\s''\ < ^ - \s^\] (4.28) 



= {s G m2 : < (j's^ < 7^, = - \s^\), \s^\ < ^ - |s^|} (4.29) 



It has eight boundary triangles given by 

where I, J = 1, 2; I ^ J; a, = =b. 
We easily compute 

Volo(Oo) = ^ (4.30) 
while (remembering that the surfaces carry outward orientation if the edges are outgoing from v) 

Fi{t^j^^,) = \[a6ij + a'6n] (4.31) 

Labelling the edge dual to Y{t^^^,^) by ejaaa we find for the expectation value 

Volo(i;) 



'48-64 



(4.32) 

where e^'^ is the alternating symbol for /, J = 1, 2 with e^^ = 1. Expression (j4.32p is already very complicated 
to analyse for the most general edge configuration and again we refer to [18j for a comprehensive discussion. 
However, for the case of the graphs constructed in section [3] the situation becomes simple enough. Namely 
in this case the eight edges e/o-o-' have the property that ei^sigma,a' and e/^_o-,-o-' are analytic continuations 
of each other. This implies that ei^a,a'{0) = o"'e/^o-o-',+ (0) where eiaa'i^) = v is the common starting point 
of all edges. Since e(e, e' , e") = sgn(det(e(0), e'(0), e"(0))) is completely skew in e, e', e", in this case we can 
simpUfy ([02]) to 



Volo(w 



^i^h,aia[,+ -> 6/2,0-20-2.+ ' 6/3,030;;, + ) 



'48-64 

Il,/2,-f3 = 1.2 
<^l,<^2,<^3'''l''^2'°"3=* 



X [c7icj;cJ2CJ^e^i^2 + o-icr'ias^e^s/i + (T2a'^aza'./^^^]\ (4.33) 
Since (14.33P only depends on ai = cr/o-j, after proper change of summation variables, (|4.33p turns into 



48 - 8 

A,^2,/3 = l,2oi,02,03=± 



Volo(t') = Jl^^ e(e/i,oi,+,e/2,02,+,e/3,<x3,+) [aia2e^^^^ + aio^^e^-^^^ + Uia^e^^H 

(4.34) 
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Using that e(e/j^CTi,+) 6/2,^2,+ ) 6/3,0-3,+) and aia2e^^^^ are both antisymmetric under the simultaneous exchange 
(fji/i) <-> (0-2/2) etc. we may further simphfy ()4.34p to 

Yolo{v) = J\^^ Yl E ^1^2 e(ei,<xi,+ , 62,02,+, e73,,^3,+) + ^ cJ2f73 e(e/,,<^i,+,ei,^2,+, 62,03,+) 
y oi,o2,o3=± I3 ii 



,01, + 5 6/2,02,+ 5 62,03,+ )]| (4.35) 

I2 

Carrying out the respective sums over /i, /21 ^3 and using that e(e, e', e") is completely skew we can bring all 
orientation factors into one of the two standard forms e(ei,(jj,+, 61,0-2,+, 62,03,+) ™^ 6(62,01,+, 62,02,+, 61,03,+) 
respectively. After proper relabelling of the aj we find that 

Volo(i^) = Jlj^ Yj ^3 [cr2 e(ei,oi,+, 6i,<^2,+, 62,03,+) e(62,oi,+, 62,02,+, 61,03,+)]! (4.36) 
y oi,o2,o3=± 

Since e(e/,CTi,+, 6/,o-2,+, 6j,o-3,+) is skew in ai,a2 the sum over (T2 collapses to the term a2 = —ox and (|4.36p 
becomes 



Volo(v) = X] CTl [-e(6l,oi,+ ,6l,_oi,+ , 62,03,+) + e(62,oi,+ ,62,-c 



,+ , 61,03, + J 



0l,03=± 



= JIy^Y ^3 [-e(6l,+,+ ,6l,_,+ ,62,o3,+) +e(62, + ,+ ,62, + ,+ ,6l,o3, + )]| (4.37) 
V 0-3=± 

Finally, using e(6/,+,+, 6/,__+, 6j,03,+) = 0-3 e(6/,+,+, 6/,_,+, 6j,+,+) and 
e(6i,+,+,6i,_,+ ,62,-,+) = e(62,+,+,62,-,+,62,+,+) = 1 wc find 

Volo(^) = ^ (4.38) 

which does not agree with (j4.30p . 
4.3 Discussion 

Interestingly, for both valence n = 4 or n = 8 the expectation value is larger than the expected value with 
the same ratio 3/(2\/2). In general, for generic edge configurations and for higher and higher valence the 
expectation value will probably also be larger in ratio than the expected volume. This is because for a 
vertex of valence n the number of ordered triples of edges contributing to the expectation value is given 
by (g) and for appropriate choice of the orientation factors, these terms all contribute with the same sign. 
Such a choice is always possible up to topological obstructions discussed to some extent in [18]. For large n 
the polyhedron dual to the vertex will approach more and more a sphere triangulated into n polygonal faces 
of typical unit area At: /n. Hence we expect the leading n behavior of the expectation value to be given by 

^J I 71,3/6 (47r/n)3 = ^Js^^^JQ = 47r/3y^37r/4 while the expected volume should approach 47r/3. 

Surely, we have not shown that for graph topologies different from a cubical one the expectation value 
of the volume operator with respect to the dual cell complex coherent states cannot be matched with 
the classical volume value. This is because one can allow degenerate triples which decrease the volume 
expectation value. However, the discussion reveals that the question for which graphs the expectation value 
comes out correctly is far from trivial and even for natural choices the only admissible graph topology is the 
cubical one. 

Notice that the expectation value is insensitive to the embedding of the graph relative to the dual cell 
complex as long as the graph is dual to it. For non dual embeddings or graph topologies which do not 
match the cell complex topology at all, the expectation value will be completely off the correct value. This 
demonstrates that the cut - off graph must lie within a certain class which is adapted to the cell complex. 
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5 Summary and Conclusions 



Together with the analysis in [M] we have shown that the only known states of LQG which are semiclassical 
for the volume operator must be based on cubic cut - off graphs. This looks surprising at first but can 
maybe be understood intuitively as follows: 

The volume operator is a derived operator and arises from the known representation of the flux operator 
on the Hilbert space. The derivation involves a regularization step which involves cubes surrounding the 
verticies of the graph in question on whose faces the fluxes are located. In order to take the limit in which 
the cubes shrink to the verticies and in order to make the result independent of the relative orientation 
between cubes and graphs, an averaging procedure must be applied. Hence one might be tempted to say 
that the fact that cubical topology is singled out rests on the cubical regularization. 

However, this is not the case. Namely, cylindrical consistency and background independence alone 
already fix the cylindrical projections of the volume operator up to a global constant as proved explicitly 
in [HI [9]. The constant depends on the averaging procedure chosen and on whether one uses tetrahedra 
rather than cubes in the regularization. However, consistency between volume and flux quantization fixes 
that factor }10j and rules out the operator |8]. That is to say, there is no freedom left in defining the volume 
operator and therefore the details of the regularization do not matter, it is a regularization independent 
result. Hence, the preference for cubic graphs in the semiclassical analysis must have a different origin. 

To see what it is, notice that the volume operator at a vertex involves a sum over ordered triples of edges 
adjacent to the vertex of which only the those with linearly independent tangents contribute. If the vertex 
has valence n then typically there are contributions [18]. They all contribute with equal weight (up to 
sign) which is the unique factor determined in [TO]. That constant is such that each triple contributes as if 
(the tangents of) a triple of edges spans a corresponding parallelepiped. However, it is clear that generally 
far less than (g) parallelepiped are sufficient to triangulate a (dual) neighborhood of the vertex and thus it 
is not surprising that large valence cut - off graphs will not give rise to good semiclassical states. On the 
other hand, unless the graph is cubic, even at low n = 4 the parallelepiped volume contribution per triple is 
too high for the triangulation of a tetrahedron. We have seen both effects at work in the previous section. 

This result has two implications: Either one is able to find new types of states which are not constructed 
by the complexifier method or by different complexifiers than the ones employed so far such that the correct 
semiclassical behavior is recovered also for graphs of different than cubic topology. Or, if that turns out to 
be impossible, one should accept this result and conclude that, in order that the boundary Hilbert space 
of spin foam models has a semiclassical sector, one should generalize them to more general than simplicial 
triangulations of the four manifold as advocated in |32^ [33] . 
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